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ABSTRACT 

We present a detailed analysis of the magnetic field structure found in a set 
of four general relativistic 3D MHD simulations of accreting tori in the Kerr 
metric with different black hole spins. Among the properties analyzed are the 
field strength as a function of position and black hole spin, the shapes of field 
lines, the degree to which they connect different regions, and their degree of 
tangling. Strong magnetic field is found toward small radii, and field strength 
increases with black hole spin. In the main disk body, inner torus, and corona 
the field is primarily toroidal. Most field lines passing through a given radius in 
these regions wander through a narrow radial range, suggesting an overall tightly- 
wound spiral structure. In the main disk body and inner torus sharp field line 
bends on small spatial scales are superimposed on the spirals, but the field lines 
are much smoother in the corona. The magnetic field in the plunging region is 
also comparatively smooth, being stretched out radially by the infalling gas. The 
magnetic field in the axial funnel resembles a split monopole, but with evidence 
of frame dragging of the field lines near the poles of the black hole. 
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We investigate prior speculations about the structure of the magnetic fields 
and discuss how frequently certain configurations are seen in the simulations. For 
example, coronal loops are very rare and field lines connecting high latitudes on 
the event horizon to the disk are not found at all. Almost the entire system is 
matter-dominated; the only force-free regions are in the axial funnel. We also 
analyze the distribution of current density, with a view toward identifying possible 
locations for magnetic energy dissipation. Regions of high current density are 
concentrated toward the inner torus and plunging region. Dissipation inside 
the marginally stable orbit may provide a new source of energy for radiation, 
supplementing the dissipation associated with torques in the stably-orbiting disk 
body. 

Subject headings: Black holes - magnetohydrodynamics - instabilities - stars: accretion 

1. Introduction 

Accretion onto black holes is potentially the most efficient form of energy generation 
possible; it is expected that as much as several tens of percent of the rest-mass energy 
of the accreted matter can become available for radiation (Novikov & Thorne 1973; but 
see also Gammie 1999 and Agol & Krolik 2000). The rate and character of accretion are 
regulated primarily by angular momentum transport, making the nature of accretion torques 
a fundamental problem of astrophysics. These torques result from magnetohydrodynamic 
(MHD) turbulence, driven by the magneto-rotational instability (MRl) (see the review by 
Balbus & Hawley 1998). 

Large-scale numerical simulation is the best tool wc have for exploring turbulent systems. 
The first attempts to describe the global structure of accretion onto black holes employed 
codes solving the equations of pseudo-Newtonian physics, i.e. Newtonian dynamics but with 
a gravitational potential of the Paczynski-Wiita form (Hawley 2000; Hawley & Krohk 2001, 
2002; Armitage et al. 2001; Armitage & Reynolds 2003; Machida & Matsumoto 2003). For 
obvious reasons, these could provide only an approximate description of the dynamics, and 
could not explore the effects of black hole rotation. 

With the construction of a fully general-relativistic three dimensional MHD code (De 
ViUiers & Hawley 2003, hereafter DH03), it is now possible to achieve a physical description 
of black hole accretion that is much more complete. Using this code we have simulated 
accretion flows onto black holes with four different spin parameters, a/M = 0, 0.5, 0.9, and 
0.998. These runs were labeled KDO, KDI, KDP, and KDE in the first paper of this series 
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(De Villiers, Hawley, & Krolik 2003, hereafter Paper I), which presented both a detailed 
description of how these simulations were done and an introduction to the several distinct 
physical regimes found within the accretion flows. In this paper we continue our report on 
these simulations with a discussion of the structure of the magnetic field. 

Before presenting this report, we briefly review the nature and conduct of the simula- 
tions. We adopt gcomctrodynamic units G = c = 1 (Misner, Thorne, & Wheeler, 1973), and 
express evolution times in terms of black hole mass M. We also follow the convention that 
Greek indices represent spacctimc vector and tensor components, whereas Roman indices 
represent pm^cly spatial components. The simulations analyzed in this paper were evolved in 
Boyer-Lindquist coordinates {t, r, (p, 6) on a 192 x 64 x 192 spatial grid {r x (px 6). This grid 
is described as "high-resolution" in Paper I. The radial grid runs from a point just outside 
the event horizon (the event horizon and the inner boundary vary with a/M) to r — 120 M 
using a hyperbolic cosine distribution to concentrate zones near the inner boundary. The 
polar angle 9 runs from 0.045 tt to 0.955 tt with reflecting boundaries at the polar axis. Only 
the quarter plane from < < 7r/2 was included, with periodic boundary conditions in 0. 
The temporal step size. At, is determined by the extremal light crossing time for a zone on 
the spatial grid, and remains constant for the entire evolution, as described in DH03. Step 
sizes are on the order of 10^^ Af for each of the four models discussed here. 

The state of the rclativistic fluid at each point in the spacetime is described by its density, 
p, speciflc internal energy, e, 4-velocity, [/^, and isotropic pressure, P. The rclativistic 
enthalpy is /i = 1 + e + P/p. The pressure is related to p and e through the equation of state 
of an ideal gas, P = pe (F — 1), where V is the adiabatic exponent. For these simulations we 
take F = 5/3. The magnetic field is described by two sets of variables. One is the constrained 
transport (CT) magnetic field, = [ijk\Fjk, where [ijk\ is the completely anti-symmetric 
symbol and Fj^ are the spatial components of the electromagnetic field strength tensor. We 
also define the magnetic field 4- vector a/Itt^'' = *F^^Ui,, where *F^^^ is the dual of the 
electromagnetic field tensor. Because the CT fleld is guaranteed to be divergence-free, it is 
the magnetic field description most readily identifled with fleld lines. On the other hand, 
the magnetic field four-vector is the one most readily interpreted as a magnetic field for 
dynamical purposes because it arises naturally in algebraically simple forms for the stress 
tensor, T^^" . The magnetic field scalar is = h^h^, and appears in the definition of the 
total four momentum, S^j, = {ph + , where W is the boost factor. Magnetic 

pressure is given by Pmag — |||^||^- We also define auxiliary density and energy functions 
D — pW and E — D e, and transport velocity — W/U^. 

For each of the four simulations, the initial conditions were basically the same: a torus 
with a near-Keplerian initial angular momentum distribution and pressure maximum at 
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r — 25 M, containing an initial magnetic field consisting of weak poloidal field loops lying 
along isodensity contours inside the torus. Each simulation was run to time 8100 M, which 
is approximately 10 orbits at the initial pressure maximum. The radial boundaries permit 
outflow only, so the simulations have only a finite reservoir of matter to accrete. However, 
even in the case with the greatest integrated accretion {a/M = 0), only 14% of the initial 
mass was accreted in the course of the simulation. After initial transients, the regions of the 
accretion flow inside the initial pressure maximum evolve into a moderately thick, nearly- 
Keplerian, highly-turbulent disk. This structure gives the appearance of rough statistical 
time-stationarity but also evolves on long timescales. 

As outlined in Paper 1, the quasi-steady state system is usefully divided into five regions: 
the main body of the disk, the coronal envelope, the inner torus and plunging region, the 
funnel wall jet, and the evacuated axial funnel. In this paper we make a more extensive 
examination of the properties and structure of the magnetic fields in these regions. Note 
that outside r ~ 25M (the location of the initial pressure maximum), we arc not simulating 
an "accretion disk" at all, as material in that region moves outward as it absorbs the angiilar 
momentum taken from matter at smaller radii. For this reason, our analysis concentrates 
on the region r <^ 25M, and especially on behavior in the vicinity of the innermost stable 
circular orbit. In §2 we examine the magnetic field strength and distribution throughout the 
flow. In §3 we quantify properties of the fleld geometry, and in §4 we determine where fleld 
dissipation would likely be important. Section 5 summarizes our findings. 



2. Distribution of Magnetic Pressure 

The late-time structure of the magnetic scalar (whose physical identification is 

relatively simple in the context of the fiuid frame-it is twice the magnetic energy density 
there) and its distribution in the accretion flow are found in Figure 1, which shows color 
contours of the azimuthal-average, at t = 8080 M for aU four simulations. The 

magnetic pressure is expressed in units of the initial maximum gas pressure in each simulation 
to permit comparison across models. (The initial gas pressure maximum increases with 
increasing black hole spin.) Gas pressure contours (at t — 8080 M) arc shown as white lines 
in order to locate the inner torus (left panels) and the main disk body (right panels). The 
gradual thickening of the inner torus with increasing black hole spin is highlighted by the 
gas pressure contours. The left panels show that the magnetic pressure is greatest on the top 
and bottom surfaces of the inner torus, and that the regions of greatest magnetic pressure 
overlie regions where the gas pressure abruptly drops. In the funnel region, the magnetic and 
gas pressure proflles are very nearly circular, though magnetic pressure is dominant in this 
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region. Although the magnetic pressure distribution looks irregular, the total pressure, gas 
plus magnetic, is relatively smooth; total pressure is predominantly cylindrical in distribution 
outside the funnel and spherical within (see fig. 6 of Paper I). Figure 1 shows that, near the 
black hole, magnetic pressure tends to increase with black hole spin. Near r = 5 M. the 
typical magnetic pressure grows by roughly an order of magnitude from a/M = to 0.998; 
near the inner radial boundary, the magnetic pressure increases by close to two orders of 
magnitude from zero to maximal spin. This stands in contrast to the magnetic pressure 
in the main disk body at large radii (r > 30 M), which does not vary much with spin. 
Figure 1 also shows that magnetic pressure diminishes rapidly both with height above the 
disk surface and with increasing radius. Gradients in magnetic pressure are very large: in 
just a few gravitational radii outside the inner boundary, the pressure falls by 1-2 orders of 
magnitude. The magnetic pressure averaged over radial shells falls roughly as from the 
inner radial boundary to the outer. 

The spatial relation between gas and magnetic pressures is shown in Figure 2, which 
presents the ratio of azimuthally-averaged gas pressure to magnetic pressure, P = 
(see also fig. 8 in Paper 1). The main disk body and inner torus show up quite clearly in 
this figure since f] is greater than unity everywhere in these regions. In the dark red regions 
of the disk body /3 can be much higher, as high as ~ 1000 in some locations, although it is 
generally lower, with a mass- weighted shell- average value of /3 ~ 10-100 at all radii within 
the main disk body. In the inner torus and plunging region, (3 tends to smaller values with 
decreasing radius for the a/M — 0, 0.5, and 0.9 simulations, but remains elevated in the 
a/M — 0.998 simulation, whose plunging region is both very compressed radially in Boyer- 
Lindquist coordinates and not well-resolved in our simulation. For all but the a/M = 0.998 
simulation, the mass-weighted mean that is 10-100 in the disk body drops to ~ 2-6 near 
the inner radial boundary, while the volume-weighted mean f3 falls from ~ 0(10) near Vms 
to 0.3-1 at the inner boundary. This decrease is not evenly distributed, however. Through 
the inner torus to the plunging region, the density contours focus toward the equatorial 
plane and contours of elevated /3 tend to follow. Although the aspect ratio (H/r) of gas 
pressure-dominated material narrows in the plunging region, in the equatorial plane itself 
/3 remains larger than unity all the way to the inner radial boundary. Figure 1 showed a 
relationship between magnetic pressure near the black hole and the spin of the black hole. 
Figure 2 shows that there is no comparable trend in f3. Rather, as the black hole spins faster, 
the gas pressure and density in the inner disk become larger relative to the mass stored in 
the outer disk, and the magnetic pressure grows proportionately. 

In the coronal envelope /? is of order unity; the funnel wall jet and other outflows within 
the corona have j3 > 1. Generally, (3 decreases with increasing distance away from the 
midplane. More than two density scale-heights above the plane, ~ 0.1 in most places. 
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Because of field buoyancy one expects the corona to be significantly magnetized. The local 
stratified simulations of Miller & Stone (2000) found j3 ~ 10~^ in the corona, but in their 
simulation there was relatively little mass outflow, due to the isothermal equation of state 
and the local disk approximation. In our simulations there is considerable outflow from the 
inner, hot portion of the disk into the corona, which helps to maintain j3 at more modest 
levels. The exception, naturally enough, is the axial funnel which is magnetically dominated 
and has /3 <?C 1. The centrifugal barrier effectively prevents any significant flow of disk 
material into this region. The funnel does, however, contain regions of anomalously elevated 
/3, due to shock heating of the tenuous funnel gas, which tends to enhance the pressure of 
the gas while its density remains extremely low. The final region of interest, the funnel- wall 
jet, is a region defined by two properties: the gas is unbound {hUt < —1) and outbound 
{Sr > 0) (see Paper 1). Although the jet contours are not explicitly shown in Figure 2 (see 
figs. 3 and 10 of Paper 1), they occupy a narrow band with j3 ~ 0.3 that lies between the 
substantially lower values typical of the axial funnel and the higher values of the corona; the 
jet lies poleward of the extended high-/? radial streaks that stand out in the right panels. 
More analysis of the funnel-wall jet and the lower density outflow within the axial funnel 
will be provided in a subsequent paper in this series. 

Another measure of the relative importance of the magnetic field is given in Figure 3 
which shows (|||&P/p/i), the azimuthally-averaged point-wise ratio of magnetic pressure to 
enthalpy. When the ratio is much greater than unity, the gas's inertia is of minor importance 
to the dynamics; this state is often called "force-free" (e.g., Blandford 2002). This figure 
clearly distinguishes the axial funnel, the main disk body, and the coronal envelope. It is 
clear that nowhere in the main disk body, corona, inner torus, or plunging region is the 
force-free condition met. In the main disk body and most of the inner torus, the ratio is 
very small, typically ~ 10^"^ in the main disk body, rising to ~ 10^^ in the inner torus and 
plunging region. The ratio ranges from 10^^ to 1 in the outer layers of the inflow through 
the plunging region (compare the prediction in Krolik 1999 that ||&||^/p ~ 0.1-1 in the most 
relativistic portion of the plunging region). In the axial funnel, the only region where the 
force-free condition is met, the ratio ranges from ~ 100 to 1000. 

It should be noted that the GRMHD code was designed primarily for regions where the 
matter inertia is nonnegligible. The code employs a density and energy floor that allows a 
dynamic range of about seven orders of magnitude in these variables. The code also places 
a ceiling on the Lorentz factor W {WceUmg = 4). In practice these limits come into play 
only at certain times within the axial funnel, which renders the accuracy of the dynamical 
properties within the funnel somewhat uncertain. For this reason we focus mainly on the 
general features of the field strength and topology in the funnel, and give less weight to 
detailed calculations of energy or energy flux there. 
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Azimuthally-averaged representations do not convey the entire picture. In Figure 4 we 
plot the fractional rms fluctuation level in magnetic energy density relative to its azimuthal 
mean in the a/M = 0.9 simulation at t — 8080 M. The rms fluctuation is defined for a 
quantity / by (Hawley & Krohk 2001) 



where (/)^ refers to the azimuthal average at fixed R and 6 of the quantity /. Fluctuations 
in field strength can be very large, and large fluctuations clearly delineate the main disk 
body and inner torus. In the a/M — 0.9 simulation the fractional rms fluctuation level 
in magnetic energy density relative to its azimuthal mean ranges from 0.8-2.0 within the 
disk. On the other hand, both the plunging region and the corona are far more regular: 
the fractional rms fluctuations in both zones drop to typically 0.1-0.3. These characteristic 
fluctuation levels are independent of black hole spin. 

Figures 1 — 4 show clearly that the different regions of the accretion flow have quite 
different magnetic characteristics. The funnel is magnetically dominated, and nearly force- 
free. In contrast, the main disk body and inner torus are gas pressure dominated, but their 
evolution is nevertheless controlled by the Lorcntz forces in the MHD turbulence. The inner 
edge of the disk is a region where the inflow becomes increasingly magnetized. The corona 
is mixed; while on average the magnetic and gas pressures are comparable, the magnetic 
pressure has a larger effective scale height, and there are distinct regions where one or the 
other is dominant. Nowhere, however, does the corona approach the force-free limit. 



Magnetic fields have directions, of course, in addition to magnitudes, and the directions 
may be visualized in terms of field lines. The overall geometry of the magnetic field is 
important, since field lines can link together different regions of the flow in a way that is 
not possible by hydrodynamics alone. A convenient way to summarize the possibilities is in 
terms of a modiflcation of the held line categorization introduced by Blandford (2002), who 
listed seven generic varieties of field lines in the vicinity of an accreting black hole. We adopt 
this list, but modify some of the definitions and include two additional types (see Fig. 5), as 
follows: 




(1) 



3. Field Geometry 



1. tangled field internal to main disk, 

2. loops connecting different parts of the disk through corona. 
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3. open field lines emerging from the main disk, 

4. field lines connecting the gas inside the plunging region to the main disk, 

5. high latitude field loops connecting the black hole to the main disk, 

6. open field lines from the black hole out through the funnel, 

7. open field lines from the plunging region to an outflow, 

8. largely toroidal field fiUing the corona, 

9. field lines connecting the black hole to the disk through the plunging region. 

We expect that the presence or absence of some of these field line types will be largely 
independent of the specific assumptions of a given simulation. For example, twisted field 
loops within the disk body (type 1) are the generic result of the MRI. Similarly, if there 
is accretion from the disk body into the plunging region, it is very hard to avoid creating 
field lines finking those two regions (type 4); the only issue is how far they extend and the 
field intensity associated with them. Toroidal coronal field (type 8) also appears hard to 
avoid whenever field is ejected into the corona. On the other hand, field line configurations 
such as type 2 may be more problematic. Such extended field loops have been suggested as 
a possible origin of high-energy emission from accretion disks. By analogy with events in 
the Solar corona, it is thought that the footpoints of magnetic field lines in the corona may 
be twisted in opposite directions by orbital shear, leading to reconnection (e.g., Blandford 
2002) and a flare. But do such held conflgurations arise naturally in a turbulent disk? 

Field line configurations such as types 5, 6, and 9 are of particular interest as they offer 
the possibility of extracting rotational energy from the black hole. Whether or not there are 
field lines connecting the accretion flow or the black hole's event horizon to inflnity bears 
directly on the issue of jet formation. The answers could depend on speciflc conditions of 
the accretion flow. For example, in these simulations we supposed that all fleld lines in the 
initial state were closed within the initial torus of matter. Thus, any fleld lines linking the 
disk to infinity must be created as a consequence of an outflowing wind or jet. The existence 
of winds and jets, and where they originate may, however, be very sensitive to the initial 
conditions and the detailed thermodynamics of the accretion flow. Another possibility is 
that accretion fiows drag in field lines from large radius, thereby maintaining a connection 
between the accreting matter and infinity. Whether flows of this type exist is a matter 
of conjecture and debate (e.g., Bisnovatyi-Kogan & Ruzmaikin 1976; Thorne et al. 1986; 
Lubow et al. 1994). 
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In this section we examine the overall geometry and global connectivity of the magnetic 
field in these simulations. We begin by displaying the general shape and character of the field 
lines. Next wc attempt to describe quantitatively the degree to which field links different 
regions of the flow, as well as the tangling of the fleld at small scales. We organize the 
discussion according to the five fiow regimes identified in Paper I and the generic field 
topologies given in Figure 5. 



3.1. Field Line Shapes 

To generate field line plots, we compute integral curves for which the tangent vector 
at each point on the curve is parallel to the local magnetic field. In many circumstances it 
might be desirable to describe the field lines in the fluid frame (using the 6^ code variables) 
because in some respects that is the most physical way to view them (see DH03 for a detailed 
discussion of the fluid-frame variables). However, the most direct way to trace these lines 
is to use the CT- variables, B\ in the Boyer-Lindquist coordinates (i.e., work directly in 
the coordinate frame). The CT variables are guaranteed to be divergence-free, and have 
a direct relationship to the electromagnetic field strength tensor (DH03). In view of these 
considerations, we define the local tangent vector, 

f^^B, (2) 

where are position vector components on the field line in the coordinate frame and A is 
some parametrization of the curve. Taking the norm of the magnetic field, 

\\Bf = ^ijB'B^ (3) 

where is the spatial sub-metric, we can rewrite this expression using (2) to obtain 

m? = 7.-^^ (4) 

This allows us to normalize (2) to the local proper length of the vector field, 

dx' B' 



ds \\B\ 



(5) 



where ds — dx^ dxK This set of first-order differential equations is integrated with a 
Runge-Kutta-Gill method using a discrete step length [ds — > As) that is kept fixed (As = 
0.01 M). The integration for each field line is begun at a starting point {ro,Oo,(t>o) and 
continued for a distance s' = Tq to either side of the starting point, where s' is the nominal 



-10- 



flat-space length, i.e., {ds'Y = dx'^rjijdx^ . In principle, a fourth-order Runge-Kutta scheme 
has an error that scales as (As)^. However, in this case, because evaluations of the right hand 
side of the differential equation depend on interpolations in gridded data, the error scales as 
(As)^. The fractional accumulated error in the field line's coordinates after an integration 
spanning s' = Tq should then be ~ 10~^(ro/M)~^/^(//ro)^^, where / is the physical bending 
scale. To keep this error to less than one grid cell requires that I/to ^ O.l(ro/M)~^, hence 
the conservative choice in step size As. To visualize these field lines, we take the output 
of the integrator, which consists of a fist of grid-space coordinates for each field line, and 
render these using a coordinate mapping from the Boyer-Lindquist (r, 6, 0) triple onto a 3D 
visualization grid. 

We begin our description of the field geometry with the main body of the disk, which 
consists, as anticipated, of turbulent, tangled toroidal fields (type 1). In the main disk body, 
the field geometry is controlled by two effects: orbital shear and turbulence. The former 
continually draws out radial field into azimuthal, while the latter twists the field in all 
directions on scales smaller than the local disk thickness. The resulting field line shapes are 
clearly shown in Figure 6. For the most part, the field lines wrap around the disk azimuthally 
as tightly-wound spirals, but there are numerous small twists and tangles. The magnetic 
field geometry is qualitatively the same for all black hole spins. 

In the corona, the field is generally much more regular, and at times is almost purely 
azimuthal. Field-line category 8 prevails here. There are almost no tangles; adjacent field 
lines run very nearly parallel to one another. We must qualify this statement by acknowl- 
edging that we can examine only selected timesteps from a full simulation and that there 
is variation with time. For example, there are some sampled time steps in which the field 
in the corona more nearly resembles the field in the main disk body as shown in Figure 6, 
presumably due to the passage of magnetic pressure bubbles through the corona, a feature 
clearly seen in animations of the density variable, p. More generally, however, poloidal field 
loops (type 2) appear only very rarely in the corona. When they do, it is near the surface 
of the disk, where small scale turbulent field structures emerge from the disk body. In the 
disk body, the MRl maintains turbulence, and the energy density of turbulent fluid motions 
is comparable to or greater than the energy density in magnetic field. In the corona, on the 
other hand, the energy density of the magnetic field is generally greater than the energy den- 
sity in random fiuid motions, so the field adopts a shape primarily governed by equilibrating 
magnetic forces, subject to the constraint imposed by the large-scale orbital shear. 

The inner part of the disk consists of the inner torus (where gas pressure reaches a 
local maximum) and the plunging region. The plunging region lies inside the turbulence 
edge (Krohk & Hawley 2002), that is, the point where the net infiow speed becomes greater 
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than the amphtude of the turbulent motions rather than the other way around. In the 
plunging region, the field structure is intermediate in character between the disk body and the 
corona. Because inflow starts to dominate over turbulent fluctuations, the field is primarily 
controlled by relatively smooth stretching of field lines. Frame-dragging also helps in this 
regard. However, complete inflow dominance is not achieved until the flow reaches deep into 
the plunging region. In the a/M = 0.9 simulation, the innermost stable circular orbit is at 
2.32 M, at least a small amount of tangling persists throughout the plunging region, while 
in the a/M — 0.998 simulation (not shown), where the plunging region is poorly resolved 
(Paper I), extensive tangling remains down to the inner boundary. Field- line types 4 and 9 
predominate in this region, as the plunging region is linked magnetically to both the main 
body of the disk and the inner boundary of the simulation. 

The magnetic field in the axial funnel is created in the initial accretion event when 
relatively strong field lines first reach the horizon. Material drains off the field lines and 
they rapidly expand out from the equator to fill the funnel. Once established, these field 
lines remain largely unchanged through the remainder of the simulation, and in polar force 
balance with the corona. These are field lines of type 6, a split-monopole configuration 
with loosely-wound helices. The dominant component is radial, but there is an azimuthal 
component that increases with black hole spin and proximity to the black hole. The specific 
angular momentum of material in this region is 2-3 orders of magnitude smaller than in the 
nearby disk, so nearly all the rotation is due to frame-dragging. Conversely, at increasing 
distances from the hole, the funnel field lines become almost purely radial. Because the 
angular momentum barrier prevents any matter located initially in the torus from finding 
its way into the axial funnel, there are almost no field line connections linking any part of 
the disk to the outflow. There are likewise no field lines that start at high-latitude points on 
the inner radial boundary, pass through the axial funnel, and lead to the disk body. That 
is, field lines of types 5 and 7 are entirely absent. 

The field lines in the funnel-wall jets (not shown) resemble those of the axial funnel, 
except that they are more tightly wound, considerably more so for the high-spin simulations, 
a/M = 0.9 and 0.998. This is consistent with the observations made in Paper I that the 
gas in the funnel- wall jets has specific angular momentum comparable to the value at the 
injection points, the surface of the inner torus. 

The effect of frame dragging on the field lines is especially prominent at the smallest 
radii, just outside the event horizon. Figure 7 makes this point dramatically. When a/M = 0, 
field lines threading the event horizon outside the equatorial plane are almost perfectly radial. 
Those passing through the horizon in the plane are swept back azimuthally because they 
are being carried in with matter having substantial angular momentum. However, as a/M 
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increases, the tightness of the azimuthal winding of all field lines, whether associated with 
accreting matter or not, grows sharply more dramatic as frame-dragging enforces rotation 
in the Boyer-Lindquist frame. 

Close inspection of the a/M — panel in Figure 7 reveals a curious effect, specific to 
that case. As the inner boundary is approached, the lapse causes the transport speeds to 
fall, reversing the sense of shear. Because the stream lines in the equatorial plane arc more 
tightly-wound than the field lines, the shear controls the correlation between and B'^, and 
the field lines bend backward. 



3.2. Field Line Connectivity 

An important property of magnetic field lines is their ability to link distant regions 
through magnetic tension. Indeed, it is the ability of the field to connect orbiting fluid 

elements that is the cause of the MRI and the resulting angular momentum transport in 
accretion flows. In the context of Kerr black holes, magnetic fields offer the additional pos- 
sibility of linking the rotating spacetime near the hole with more distant regions, potentially 
extracting rotational energy from the hole itself. 

In an attempt to quantify the length scales connected by the field lines, we define the 
wandering index, a simple measure of the range of coordinate values spanned by a given field 
line. Specifically, we define the wandering index AX for a poloidal coordinate X = {r/vo, 9) 
as the difference between the maximum and minimum X-value along a particular field line 
of length s' — To- Next, to compute mean values of this index we choose approximately 100 
fleld line starting points within a specifled range of angular coordinates on the spherical shell 
with radius Vo at each of ten different times spanning the last 720 M of the simulation. The 
probability density for the starting points is the density of fleld-lines piercing that radial 
shell; i.e., it is proportional to {^Jg^ / — g)B^ . In these calculations, we do not take account 
of fleld lines that reach the inner boundary. 

The wandering index for different coordinates takes on characteristic values in certain 
limits. For example, when the fleld line is purely toroidal, both Ar/vo and are identically 
zero. Purely radial fleld has Ar/vo = 1 and = 0. Field lines conflned to the equatorial 
plane have A^^ = 0, and the size of their Ar/vo index indicates the tightness of their spiral 
winding: Ar/vo — > is the tight-winding limit, and Arjvo = 1 is the purely radial limit. 
We caution, however, that the scaling of AX with integration path-length s' is difficult to 
estimate. 

In the disk body, a typical fleld line connects regions separated by Ar/vo — 0.2, with 
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little dependence on a/M (Figure 8). Thus, on this scale of field line length, the field 
lines are moderately wound (as seen in Figure 6). For the a/M = and 0.5 models, the 
radial distance spanned increases inside the marginally stable orbit as the flow changes from 
one that is primarily turbulent to one that is primarily spiral inflow. For the high-spin 
models {a/M = 0.9 and 0.998), inside the static limit {r static = 2M in the equatorial plane, 
independent of spin) frame-dragging ensures that most of the field line's physical length is 
stretched in azimuth rather than radius, so that Ar/ro falls. 

The wandering index for the polar direction is also shown in Figure 8. In the disk body, 
field lines move typically 0.08-0.15 radians in polar angle in the course of traversing one 
radius in length. By contrast, the typical aspect ratio of these disks as found from the gas 
density scale-height is H/r ~ 0.2 (Paper I). Thus, although the magnetic field scale-height 
is roughly double that of the gas, the vertical distance across which an individual field line 
moves in a length of one radius is less than a single gas scale-height. The focusing of the 
infiow toward the equatorial plane further reduces the polar wandering index as r decreases. 



The preceding measures mostly refiect the larger-scale properties of the field line ge- 
ometry. While a field line may not extend over significant distances, and hence may have 
a small wandering index, it may nevertheless be far from smooth and uniform. To quan- 
tify the small-scale structure of the field lines, we define a field line tangling index, ^, by 
accumulating the angle between adjacent tiny segments along field lines of length s' — To'- 



where n*^^-j is the unit three-vector in the direction from point k — 1 to point k along the 
field line and r/jj is the flat space metric. For this index, the sequence of points is not the set 
of points generated at each integration step but rather a set along the integrated field line 
separated by As' = 0.025 r. We chose this scheme because the radial coordinate grids of the 
simulations were very nearly logarithmic, with separations varying between Ar = 0.025 r for 
the a/M — case to Ar = 0.028 r for the a/M = 0.998 run. By choosing this sampling 
scheme, we assure ourselves of achieving a fixed resolution level across the different regions 
of the simulations. 



3.3. Field Line Tangling 




(6) 



fc 



In the limit that the resolution scale is small compared to the length scale of large- angle 
field line bending, the tangling index has a simple interpretation. For small individual angles. 
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(6) becomes 




1/2 



(7) 



Vij '^(fc+1) ^(fc) 



k 

When the individual bend angles are small, n*^^-) is effectively a continuous function of distance 
along the field line and so may be expanded as a Taylor function in s'. Taken to second-order, 



That is, the tangling index is the ratio of the length of the field line to the harmonic mean 
of its radius of curvature. In the unlikely limit of a perfectly straight field line, ^ = 0. For a 
circular segment in any plane, ^ = 1 (including the equatorial plane, where such a segment 
would be purely toroidal) . Because we compute a discrete approximation to the integral with 
40 (i.e. 1/0.025) individual angles, the maximum possible tanghng index is ^ 407r ~ 125. 

The results of computing this measure arc displayed in Figures 9 and 10. Over one radius 
of field line length, the accumulated small-scale bends in the disk body sum to 4 - 7 radians, 
with a tendency for the tangling to increase with increasing a/M. In approximate terms, the 
distribution of individual bend-angles dN/da inside the disk is proportional to a"^, so the 
contribution to the total bending from each decade of angular scale is roughly constant, and 
the cumulative bending distribution rises logarithmically with bend angle (Fig. 10). Such 
a result might have been expected given the turbulence in the disk. Inside the marginally 
stable orbit, the degree of small-scale tangling drops sharply as turbulence gives way to a 
comparatively smooth flow. The deeper inside the plunging region the flow reaches, the less 
tangled the field lines. In the coronal region the field lines become much smoother, with 
the tangling index falling to 1-2. In the limit of a perfectly smooth azimuthal field line, the 
index would be identically 1 as described above; thus, the coronal field lines are very smooth 
indeed. 



(8) 



where R{s',j^^) is the local radius of curvature. Thus, in the limit As' 0, 




(9) 



4. Other Field Properties 



4.1. Symmetry of B'^ 



Since the initial condition consists of purely poloidal field loops centered on the equator, 
the initial radial field is antisymmetric across the equator, as is the sign of B'^ that results 
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from subsequent shear. This holds true through much of the evolution. During the first two 
orbits the total toroidal flux grows rapidly and antisymmetrically. However, the accretion 
flow that forms is unstable to vertical oscillations across the equator, breaking the perfect 
antisymmetry. As a result, the net toroidal flux in each hemisphere gradually declines as 
matter initially in one hemisphere is mixed into the other. The global net toroidal flux can 
depart from zero once field begins to enter the black hole, but it never becomes as large as 
the net flux in one or the other hemisphere. At the end of the simulations there are still some 
obvious antisymmetries in B'^ in the corona, and near the funnel in the a/M — 0.9 and 0.998 
simulations, but no systematic antisymmetry remains within the disk itself. At the end of 
the simulations, the net B'^' flux compared to the total, i.e., J B'^^/^d^x^/ J \B'^\^J^ dPx^ is 
~ 1% when a/M = and 0.998, 10% when a/M = 0.5, and 8% when a/M = 0.9. 



4.2. Correlation with velocity 

In computing idealized axisymmetric, smooth, and time-steady accretion solutions with 
large-scale magnetic flelds, Li (2003a,b) shows that if the ratio {h't'u'')/{h''u't') = 1, MHD 
accretion onto black holes can occur with zero electromagnetic transport of energy. Presum- 
ably this conclusion is closely linked to the fact that when this ratio is unity, the velocity 
and magnetic flelds in the disk plane are parallel. In the ideal MHD limit, where E = vx B, 
this means E = 0. When that is so, the Poynting vector must also be identically zero, and 
no energy is transported electromagnetically in the course of accretion. We can explicitly 
measure this ratio with data from our simulations. 

As illustrated in Figure 11, in general this ratio is quite different from unity. In most of 
the corona, disk body, inner torus, and plunging region it falls in the range —0.3 — hO.3, but 
there are locations near the funnel wall where the ratio spans a considerably wider range, 
with thin regions in the vicinity of the jet having values near unity. Values near unity are 
also found at the surface of the plunging flow above and below the main equatorial inflow, 
but even there the ratio fluctuates substantially in space and in time. We therefore conclude 
there is nothing to prevent signiflcant electromagnetic energy fluxes in these accretion flows. 



4.3. Dissipation Regions 

These simulations were conducted with no explicit resistivity, and, consequently, they 
do not address where the magnetic fleld is dissipated or at what rate. However, regions of 
high current density are candidates for regions of high magnetic dissipation because high 
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current density may trigger anomalous resistivity through mechanisms such as ion-acoustic 
turbulence (e.g., as suggested for the Solar corona by Rosner et al. 1978). We can locate 
these candidate dissipation regions by computing the current density from our simulation 
data. The current 4-vector is given by 



where is the covariant derivative and F^" the electromagnetic field-strength tensor. The 
covariant derivative simplifies to a simple derivative using the anti-symmetry of F'^'^ (eq. 23 



where Fx^ is directly related to the CT magnetic field and EMFs (eqs. 14 and 35 of DH03), 
the EMFs being simple functions of the CT magnetic field and transport velocities V^. 
Evaluation of using code variables requires data for three adjacent time steps, since 
time derivatives must be evaluated. The current density then follows directly from these 
calculations, ||J|p = J^- 

Figure 12 shows azimuthally-averaged || J|p at a late time in each of the four simulations. 
The values are normahzed to the initial torus total energy, Eq (= J d^x^ ^J^T^ in the 
initial state), to facilitate comparison between simulations. In all cases, regions of high 
current density in the accretion flow are found in extended sheets that run roughly parallel 
to density contours in the inner torus and plunging regions; currents in the main disk body 
are significantly weaker and less well organized, while the coronal envelope is quiescent. The 
axial funnel shows a current distribution similar to the distribution of magnetic pressure 
(Fig. 1). In all cases, the dynamic range between || J|p in one of the sheets and in the adjacent 
inner torus material is often two orders of magnitude. As with the magnetic pressure, there 
are also several systematic trends with black hole spin. The vertical thickness of the region 
in the disk body where high current density is found increases with a/M. The absolute level 
of II J|p (even after normalization to Eq) also increases with a/M, climbing by about a factor 
of 300 from a/M = to a/M = 0.998. 

Figure 13 shows || J|p/£'o in the equatorial plane (^ = 7r/2) at a late time in each of the 
four simulations. This figure emphasizes the spiral structure of the currents. Here again, the 
steep radial gradient in current density can be seen. The regions of most intense current are 
almost all found within ~ 2rms- 



of DH03) 




1 



(11) 



If we knew how to relate current density to dissipation rate, these maps would provide 
a basis for a map of heating in the accretion flow. In the absence of a physical model for 
the relation between current density and dissipation, we can instead give a qualitative sense 
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for the heating distribution by integrating \\J\\^ over volume in the disk, i.e., computing 

J' = J a/— ^ II but excluding those regions with p < 10~^pmax, where pmax is the 
greatest density found anywhere in the simulation at that time. We ignore low-density 
regions in order to concentrate on possible heating in the accretion flow proper. The results 
are shown in Figure 14, again normalized to Eq. The left panel shows the integrated current 
density over the entire range of radii; the right panel shows the same quantity plotted as a 
function of r /vms- Both panels show the systematic increase in current density with black 
hole spin, particularly for r < lOM. In the disk body, J declines steeply with increasing 
radius. Between and lOr^s, —2 < d\n. J / d\n.r < —1, and the slope steepens at still 
larger radii. In the a/M — simulation, the integrated current density levels off slightly 
through the marginally stable orbit before rising sharply again near the inner boundary. For 
the other three simulations, the integrated current density has a sharp break upward near or 
inside r^s, rising one to two orders of magnitude as the inner radial boundary is approached. 
Since frame dragging has been shown to increase magnetic pressure and azimuthal stretching 
of the field lines, this dramatic increase is not surprising. If the current density is indicative 
of the dissipation density, much of the total dissipation, especially around the most rapidly- 
spinning black holes, could occur very deep in the relativistic potential. 



5. Discussion and Summary 

In Paper 1 we divided the accretion flow into five distinct regions, each with its own 
characteristics, as follows: the main disk body, the inner torus and plunging region, the 
coronal envelope, the funnel, and the funnel wall jet. In this paper we have examined the 
magnetic field strengths and topologies in these regions. 



5.1. Field Characteristics in the Different Flow Regions 

The main disk body is characterized by MHD turbulence due to the MRI. As has 
been shown in previous studies, the magnetic field is very tightly wound by the differential 
rotation, but with significant poloidal tangling due to turbulence. The field wandering index 
shows that the poloidal field lines extend farther in the radial direction than in the polar. 
The field energies are subthermal, with plasma /3 fa 10-100 on average. Strong fiuctuations 
are the rule, however, and /3 can exceed 1000 in some portions of the main disk body. 

In the plunging region, the fiow transitions from turbulence-dominated to spiraling 
infiow. A consequence of this is that the tight field winding loosens somewhat as the radial 
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component of the velocity grows relative to the azimuthal component. This behavior is 
quantified approximately by the field-wandering and field-tangling indices, which show (see 
Figs. 8 and 9) that there is a sharp change of behavior near the marginally stable orbit. 
Field lines inside that radius travel farther in radius and much less in polar angle than in the 
disk body. This is a consequence of the transition from turbulent motions to more regular 
spiral infiow; that is, this break is yet another way of marking the turbulence edge defined 
by Krohk & Hawley (2002). Inside the plunging region the field is amphfied relative to the 
fluid pressure. Near the horizon above and below the equator, (3 falls from ~ 10 near r^s to 
less than 1. The equatorial flow itself maintains /3 ~ 10 or greater throughout the plunging 
region. These conclusions apply most clearly to the a/M = 0, 0.5, and 0.9 simulations; the 
plunging region in the a/M — 0.998 simulation is less well-resolved. 

The fleld in the corona is, on average, in equipartition with the thermal energy there, 
and is dominated by the toroidal component; fleld line plots look smooth compared to those 
in the disk body. There are relatively extended poloidal fields, but these are more associated 
with outflowing toroidal coils. When field emerges from the disk body it carries some of the 
turbulent tangling with it, but this seems to be rapidly smoothed out. Gas is carried into the 
corona as well as magnetic field, so the corona does not become a magnetically dominated 
force-free region. The magnetic field does have a larger scale-height than the gas pressure, 
though. 

By contrast, the field in the axial funnel is essentially radial outside the ergosphere, 
with a toroidal component that becomes larger with increasing black hole spin. In the 
funnel the magnetic field energy dominates over the matter energy, and the field can be 
regarded as force-free. The magnetic energy is not large compared to the magnetic energy 
in the remainder of the flow, however. The magnetic pressure in the funnel is in equilibrium 
with the total pressure in the corona and inner disk. From the fleld-connectivity data, we 
see that the magnetic fleld structure may be divided into two almost wholly independent 
regions: the accretion flow proper and the axial funnel. Field lines passing through the 
funnel almost never enter the other regions. 

In the funnel- wall jet, the magnetic fleld conflguration resembles that of the axial funnel, 
but the fleld lines are more tightly wound, consistent with the observation that the gas in 
the jet has speciflc angular momentum comparable to that at its injection point, above the 
inner torus. The jet is also less strongly magnetized than the axial funnel, with (3 ~ 0.3. 

The categorization of fleld line types listed in Table 1 and illustrated in Figure 5 provides 
a handy summary of the global fleld structure. It is important to recognize that only a subset 
of all possible field line types are found with any frequency in these simulations; the table 
indicates which these are. One of the most important global properties of the field is its 
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inter-region connectivity. Prom this point of view, the different regions fall into two groups: 
the main disk body is linked to both the corona and the inner torus, while many field lines 
stretch from the inner torus to the plunging region. On the other hand, the field lines in the 
axial funnel have essentially no connection to the other regions. 

5.2. Astrophysical Implications 

The simulation results provide some insights into potentially important astrophysical 
processes that might occur in black hole accretion. For example, Livio et al. (1999) argued 
that energy extracted clectromagnetically from the accretion disk proper would always out- 
weigh that extracted from spinning black holes because the product of magnetic stress and 
area would be larger for the disks than for black hole event horizons. As Figure 1 shows, the 
total intensity of the magnetic field increases sharply inward. This is as true of the poloidal 
portion as the toroidal. Increases by factors ~ 30 from r = 5M to the plunging region are 
typical. Although it is true that, outside the funnel region, there is no large-scale poloidal 
field in this simulation, such large contrasts in field strength from disk to plunging region 
raise doubts about the Livio et al. hypothesis. Moreover, the force-free field in the evacu- 
ated funnel seems best configured to permit processes closely related to the Blandford-Znajek 
mechanism (Blandford & Znajek 1977), while the strength of the funnel field is comparable 
to that in the plunging region. These results suggest that accretion could well lead to a 
strong enough field near the event horizon that a genuine Blandford-Znajek process might 
yield more energy than competing processes associated more directly with the main disk 
body. 

It is often asserted (in fact, frequently in connection with studies of the Blandford-Znajek 
mechanism) that the magnetospheres of black holes should be, in large part, effectively force- 
fi-ee (e.g., Okamoto 1992, Ghosh & Abramowicz 1997, Ghosh 2000, Park 2000, Blandford 
2002, Komissarov 2002). In our simulations we find that the only region where the force- free 
approximation is applicable is the funnel. In the disk proper and most of the corona, the 
field structure is thoroughly matter-dominated; in the outer corona and parts of the plunging 
region, the energy density in the magnetic field begins to approach p/i, but never exceeds it. 
Thus, the relativistic MHD approximation may be better-suited to work in this area than 
the relativistic electrodynamic approximation. It is possible, though, that cooling would 
alter the properties of the corona somewhat by significantly reducing the gas pressure scale 
height relative to that of the magnetic field. 

A common suggestion for the origin of high-energy emission from accretion disks around 
black holes is that the footpoints of magnetic field lines in the corona are twisted in opposite 
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directions by orbital shear, leading to reconnection (e.g., Blandford 2002). This process is 
initiated by establishment of a radial magnetic field link between two fluid elements arranged 
so that the outer one is advanced in orbital phase relative to the inner one. As the inner 
fluid element catches up to and passes the outer, an X-point would be formed. We see 
little evidence of this sort of process. What appears to happen instead is that when radial 
magnetic fleld links are created, it is between two fluid elements at nearly the same azimuthal 
angle (the closer two points are, the more likely they are to be magnetically connected). The 
orbital shear then creates a long azimuthal field line as the fiuid element at the smaller 
radius moves ahead of the other fiuid element. If the inner fiuid element moves a full 27r in 
azimuth ahead of the outer one before reconnection occurs, the field lines wrap parallel to 
one another, rather than crossing near their footpoints. 

It is generally believed that regions of high current density are likely to promote anoma- 
lous resistivity, and therefore be markers for regions of concentrated heat dissipation. Al- 
though we find that, if anything, current sheets are especially rare in the corona, they are 
very common and very intense in the innermost parts of the accretion fiow, particularly the 
plunging region. The possibility of substantial heat release inside r^s has an interesting 
consequence. In the Novikov-Thorne theory of relativistic accretion disks, the energy made 
available for radiation is the difference between the potential energy lost by accreting mat- 
ter and the work done by accretion stresses. Because they assumed that all stresses vanish 
inside rms, there is no energy released inside that radius. If the Novikov-Thorne model is 
generahzed to allow non-zero stress at r^s (e.g., Agol & Krohk 2000), the amount of energy 
dissipated inside the disk rises, but it is still determined in essentially the same way. Dissi- 
pation inside rms is an entirely new and independent mechanism of energy release. It is not 
necessarily associated with angular momentum transport; it could result from entirely local 
processes (e.g., Machida & Matsumoto 2003). Nonetheless, if the time to radiate this heat 
is shorter than an infall time, the photons released could carry off the energy to infinity, 
adding to the light seen by distant observers. 

One of the main purposes of carrying out these simulations across a broad range of 
black hole spins is to extract information previously unobtainable from pseudo-Newtonian 
simulations of accreting tori. The spacetime rotation caused by spinning black holes is a 
notable example of this sort of effect. Perhaps the most striking spin-dependent property of 
the magnetic fields is the increase in magnetic pressure near the black hole with increasing 
spin. This pressure increase is found at the surface of the inner torus and plunging region (i.e., 
off the equatorial plane), as well as deep in the axial funnel. There is also a corresponding 
increase in gas pressure in this region, so that the ratio of pressures (/3) remains effectively 
unchanged. The growth of magnetic pressure is due in part to the dragging of field lines by 
the black hole, and in part to the inward displacement of the marginally stable orbit with 



- 21 - 



black hole spin, which tends to allow the MRI to operate closer to the event horizon. There 
is also a strong spin-dependent growth in the intensity of currents in the vicinity of the black 
hole. In the a/M = 0, 0.5, and 0.9 simulations, the high current regions are organized in 
sheets through the inner torus; these current sheets also have a spiral character, inherited 
from gas motion in the vicinity of the plunging region. The a/M = 0.998 simulation shows a 
much more intense and turbulent current distribution, reflecting the fact that the plunging 
region is extremely narrow in Boyer-Lindquist coordinates and the the inner torus extends 
deep into the ergosphere. This strong spin-dependence of the currents suggests that magnetic 
dissipation may be especially intense in accretion flows near extreme black holes. 
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Table 1. Field Line Varieties (after Blandford 2002) 



TY]:)e Numl)er 


Descri]:)tioii 


Common'? 


1 


tangled within the disk body 


yes 


2 


loops linking regions of disk through corona 


no 


3 


poloidal from disk to infinity 


no 


4 


linking accreting matter in plunging region to disk body 


yes 


5 


linking high latitudes on event horizon to disk body 


no 


6 


poloidal from event horizon to infinity 


yes 


7 


poloidal from matter in plunging region to infinity 


no 


8 


smooth tightly-wrapped coronal spiral 


yes 


9 


poloidal field linking horizon to accretion fiow 


yes 
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Fig. 1. — The azimuthally-averaged magnetic pressure, ), at t = 8080 M normalized 

to the maximum initial gas pressure, Pmax, in simulations with a/M = (first row from top), 
0.5 (second row), 0.9 (third row), and 0.998 (fourth row). The left column shows the region 
inside r = Sr^s, plotted on a logarithmic scale that emphasizes the increasing dominance 
of magnetic pressure with black hole spin in the innermost region of the computational 
domain. The right column shows the main body of the disk, in a color scale that emphasizes 
the dominance of gas pressure in the disk body (the upper color ranges are saturated here). 
Three azimuthally-averaged gas pressure contours are also shown, {Pgas), scaled to the initial 
pressure maximum, 0.1 P^ax (dotted line), Pmax (solid line), 10 Pmax (dashed line). 
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Fig. 2. — The azimuthally-averaged ratio of gas to magnetic pressure, (3 = (-P)/(|||&||^), at 
t = 8080 M in simulations with a/M = (first row from top), 0.5 (second row), 0.9 (third 
row), and 0.998 (fourth row). The left column shows the region inside r = Sr^s, the right 
column shows the main disk body; all are plotted on the same logarithmic scale. Three 
azimuthally-averaged density contours are also shown, (p), scaled to the initial pressure 
maximum, 0.1 pmax (solid line), 10^^ Pmax (dotted line), 10~^ pmax (dashed line). 
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Fig. 3. — The azimuthally-averaged ratio of magnetic pressure to enthalpy, {\\\h\\'^ / {p h)) , 
aX t = 8080 M for the simulation with a/M = 0.9. The left panel shows the region inside 
r = 3 Tms, the right panel shows the main disk body; both are plotted on the same logarithmic 
scale. Three azimuthally-averaged density contours are also shown, (p), scaled to the initial 
pressure maximum, 0.1 pmax (solid line), 10""^ pmax (dotted line), 10~^ pmax (dashed line). 
Images from the other three simulations are very similar. 



Fig. 4. — The fractional fluctuation in azimuth for at t = 8080 M in the simulation 
with a/M = 0.9. The left panel shows the region inside r = Sr^s, the right panel shows 
the main disk body; both are plotted on the same linear scale. Images from the other three 
simulations are very similar. 




Fig. 5. — Schematic diagram (after Blandford 2002) showing possible magnetic hnkages 
between different regions of the accretion flow and the black hole. 
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Fig. 6. — Sample field lines in four regions of the a/M = 0.9 simulation. Upper left panel 
is the disk body (i.e., 0.457r < 6 < O.SStt, = lOM); upper right is the corona (0.327r < 
9 < O.STtt, To = lOM); lower left is the plunging region (0.4757r < 9 < 0.5257r, = 2M); 
lower right is the axial funnel (O.lvr < 9 < 0.27r, Vo = 5M). All four are at t = 7760M. 
Background colors are contours of log(p), calibrated by the color-bar. 




Fig. 7. — Sample field lines threading the innermost radial boundary in simulations with 
a/M = 0., 0.5, 0.9, and 0.998. 
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Fig. 8. — Wandering index (see text for definition) as a function of radial coordinate of 
the field line's starting point for each of the four black hole spins considered. The curves 
are: a/M — simulation (solid fine); a/M — 0.5 simulation (dotted fine); a/M — 0.9 
simulation (dashed fine); and a/M — 0.998 simulation (dash-dot fine). Left-hand panel: 
Radial wandering. Right-hand panel: Polar angle wandering. 
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Fig. 9. — The field-tangling index \E' as a function of inside the disk body (0.457r < 
9 < 0.557r: left panel) and 9o (lOM < r < 15M: right panel). The curves are: a/M = 
simulation (solid fine); a/M — 0.5 simulation (dotted line); a/M — 0.9 simulation (dashed 
line); and a/M — 0.998 simulation (dash-dot line). 
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Fig. 10. — The cumulative distribution of tangling ^(< a) = da'a'dN/da' averaged over 
ten late-time snapshots in the a/M = 0.9 simulation at the radii shown and within O.OStt 
radians of the equatorial plane. Here dN/da is the number distribution of bend angles. 



Fig. 11. — The azimuthally-averag ed ratio (fe*^ u'') / {b'' u*) , at t = 8080 M in simulations with 
a/M = (first row from top), 0.5 (second row), 0.9 (third row), and 0.998 (fourth row). The 
left column shows the region inside r = Sr^s, the right column shows the main disk body; 
all are plotted on the same linear scale. Three azimuthally- averaged density contours are 
also shown, (p), scaled to the density maximum, 10~^ p^ax (dashed line), 10~^ pmax (dotted 
line), 0.1 pmax (solid line). 
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Fig. 12. — Azimuthally-averaged at t = 8080 M in simulations with a/M = (first row 
from top), 0.5 (second row), 0.9 (tliird row), and 0.998 (fourtfi row). Tlie left column shows 
the region inside r = 3 Vms, the right column shows the main disk body; all are plotted on the 
same logarithmic scale. Three azimuthally-averaged density contours are also shown, (p), 
scaled to the density maximum, 10~^ Pmax (dashed line), 10""^ pmax (dotted line), 0.1 pmax 
(solid hne). 
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Fig. 13. — II^IP in the equatorial plane {6 = 7r/2) at t = 8080 M in simulations with (from 
left to right) a/M = 0, 0.5, 0.9, and 0.998. In each case, the outermost radius shown is 3rms- 
All are plotted on the same logarithmic scale as Figure 12. 
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Fig. 14. — Shell-integrated current density, (see text for definition) plotted as a function 
of r/M (left panel) and r/rms (right panel) for the a/M = simulation (solid line), 0.5 
simulation (dotted line), 0.9 simulation (dashed line), and 0.998 simulation (dash-dot line). 
The integral over 9 includes only that region where p > Pmax, thereby excluding the 
funnel. 



